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Gamma-Ray Bursts (GRBs) are copious sources of gamma-rays whose interaction with a plane¬ 
tary atmosphere can pose a threat to complex life. Using recent determinations of their rate and 
probability of causing massive extinction, we explore what type of universes are most likely to har¬ 
bour advanced forms of life. We use cosmological N-body simulations to determine at what time 
and for what value of the cosmological constant (A) the chances of life being unaffected by cosmic 
explosions are maximised. Life survival to GRBs favours A-dominated universes. Within a ACDM 
cosmology, the likelihood of life survival to GRBs is governed by the value of A and the age of the 
Universe. We find that we seem to live in a favorable point in this parameter space that minimises 
the exposure to cosmic explosions, yet maximises the number of main sequence (hydrogen-burning) 
stars around which advanced life forms can exist. 


Why the value of the cosmological constant A is neither 
zero nor of the order of the Planck density (the Planck 
mass to the fourth power, Mp,) remains one of the deep¬ 
est mysteries of nature. Refs mi have argued that in 
order for observers (and thus cosmic structure) to exist, 
the value of A could not be larger than 10 _12O Mpj. This 
was the first indication that the value of A could not be 
arbitrary and that requiring the existence of observers 
bounded A from above. Subsequent works (e.g., mi) 
have firmed up this argument. However, to this date no 
argument has been given to provide a lower bound to 
A; in particular, it is not clear why A does not simply 
vanisf0 Interestingly, the necessity to avoid massive life 
extinction events by GRBs can shed a new light on this 
issue. 

GRBs are potentially catastrophic events for biological 
organism^ In particular, copious flux of 7-ray photons 
with energies above 10— 100 keV could destroy the ozone 
layer of a - habitable, Earth-like- planet, exposing living 
organisms to damaging UV radiation and compromising 
its habitability. This has led to the suggestion ISCDP that 
galactic GRBs have been responsible for some mass ex¬ 
tinction events on Earth. Yet the rate or energy of nearby 
GRBs were not sufficient to avoid the emergence of ob¬ 
servers. However, such GRBs take place more frequently 
at the inner parts of the Milky Way and may cause a 
serious problem for development of life there [12j . On 
Earth and in general, in the outskirts of large galaxies, 
the most luminous GRBs, the ones around the knee of 


1 A lower bound to A based on the stability of atoms [8j, applies 
only in the distant future for spatially flat universes. 

2 In what follows when we refer to life we consider biological or¬ 
ganisms which are sufficiently complex to act as observers. Our 
considerations will be for Earth-like planets where the UV pro¬ 
tection provided by the atmosphere is due to an Ozone layer. 


the luminosity function see Fig. 3 in [T2], pose the great¬ 
est threat for the development of complex organisms [^is 
they could cause catastrophic damage even if located in 
a sufficiently nearby satellite galaxy. 

The rate of GRBs within a given galaxy depends on the 
metallicity: most GRBs take place where metallicity < 
0.3 solar and the stellar mass of the host galaxy is above 
5 x 10 7 M 0 [ill, DJ- Such low-metallicity environments 
are rare within the Milky Way and the ~ 50kpc region 
around it. Small mass, low metallicity, Magellanic Cloud 
(SMC and LMC)-type galaxies are the typical host of 
GRBs and thus the most likely location for potentially 
damaging nearby GRBs. We explore next the rate of 
such catastrophic extra-galactic events. 

The observed global GRB rate is 10 9 GRB/Gpc 3 /Gyr 
[16j^) GRB hosts have stellar masses between 5 x 10' 
and 10 10 Mq P31GL5J. Integrating the stellar mass func¬ 
tion of m in this mass range we find a stellar density 
of 10 16 M e /Gpc 3 , yielding a rate of lO" 7 GRB/M 0 /Gyr. 
This rate depends only on stellar physics and thus is inde¬ 
pendent of cosmology. Integrating (out to 200 kpc from 
the center) the dark matter profile from CDM simulations 
[15] , the dark matter mass in satellites is 20% of the total 
halo mass (2 x 10 12 ) of the Milky Way. Since ~ 1% of this 
mass is in stars |TS] we obtain a stellar mass of 4 x 10 9 in 
the satellite^] Thus, we expect 400 GRB/Gyr in satel¬ 
lites. Using the observed GRB luminosity function [16] 


3 A small fraction of mass extinction from short-GRBs m won’t 
change our results, as it would just add a baseline that equally 
penalises all values of A. 

4 We use here the rate of GRBs beamed towards a given observer. 
The overall rate is larger by a factor corresponding to the beam¬ 
ing. However, this is not relevant for this work 

5 This number is in excellent agreement with direct integration of 
the observed stellar mass in Milky Way satellites ESI 
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assuming an effective duration of 10s, we expect: 280 
(70%) > 10 52 erg; 72 (18%) > 10 53 5 erg; 4(1%) > 10 54 5 
and 0.5 — 1(0.15%) with energy of ~ 10 55 erg. Ref |T2] 
has shown that most likely there has been one GRB dur¬ 
ing the last Gyr with a fluence on Earth of 100 kJ/m 2 ; 
this fluence is the value found by El M for massive 
life extinction to take place. This event is believed to 
have caused the Orvidician extinction[23], which wiped 
out 85% of all species present on Earth at the time. Fol¬ 
lowing E1E2 we take a fluence of 100 kJ/m 2 to be the 
threshold: higher fluence would have catastrophic conse¬ 
quences for having observer^] Then the equivalent dam¬ 
aging distance from the center of the galaxy for a planet 
at 10 kpc from the center is 17kpc for 10 53 " 5 erg, 27 kpc 
for 10 54 ' 5 and 50 kpc for 10 55 erg. This implies, conser¬ 
vatively, that a region of 20 kpc (from the center of the 
host halo) should be devoid of GRB-hosting sub-halos 
for harbouring planets suitable to support observers. If 
there are no satellites in this region then there will be 
no damaging GRBs (of fluence 100/c J/m 2 ) with a rate 
much higheiQ than that on Earth. Since there is a non¬ 
zero probability of a > 10 55 erg GRB/Gyr in the satellites 
which yields a lethal 100k J/m 2 fluence at a distance of 
50kpc, we also discuss this case. 

The existence of many nearby satellites will have an 
additional effect. Many of these satellites in-fall into the 
main galaxy, bringing low metallicity material and trig¬ 
gering further star formation that will increase GRB pro¬ 
duction. We do not consider this enhanced rate in this 
work but it clearly makes development of intelligent life 
in cosmologies with numerous nearby satellites even more 
difficult than considered here. 

Any inference of cosmological parameters ought to take 
into account the selection effects which have led us to 
observe the Universe from this particular vantage point 
[25] . Of particular interest is the value of A, which can 
govern the growth of regions around large galaxies de¬ 
void of LMC-type satellites. The accelerated expansion 
induced by a cosmological constant slows the growth of 
cosmic structures, and increases the mean inter-galaxy 
separation. This reduces the number of nearby satellites 
likely to host catastrophic GRBs. Below we quantify this 
effect. 

Using iV-body simulations, we search for halos with 
dark matter masses 10 11 ' 5 < M/Mq < 10 12 that have 
no satellites of dark matter mass > 2 x 10 6 7 8 * Mq within 
a radius of 20 (50) kpc (proper); we refer to these ha¬ 
los as isolated. In particular, we use the Millennium- 


6 The amount of ozone depletion and DNA damage scales slowly 
with fluence: they are reduced by factors of 2 and 2.5 respectively 
by reducing fluence from 100 to 10 kJ/m 2 m- 

7 About 8% of the satellites are within a radius of 20kpc for a A = 0 

cosmology. Regardless of the cosmology, almost all satellites are 

within a radius of 50 — 100 kpc. 
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FIG. 1. Number of isolated halos as a function of redshift 
that are protected from the damaging effect on life of GRBs; 
we show four cosmologies and two isolated radii of satellites: 
20 kpc (upper panel) and 50 kpc (lower panel). To guide the 
eye, dashed lines show extrapolations to z = 0 due to limita¬ 
tions in re-mapping the Millenium-II simulation as this would 
have needed it to be run into the future. The number of life- 
protected regions is significantly larger in the A—dominated 
cosmologies. This is due to the effect of A at late times at 
clearing and smoothing virialised regions. Also remarkable is 
the fact that only below 2 ~ 1 the number of isolated halos 
grows significantly. This is about 7 Gyr ago, not too dissimilar 
from the age of Earth. 

II E§3 publicly available dark-matter-only simulation [^] 
through the Millennium Database portal created by the 
Virgo Consortium m■ We search for these isolated ha¬ 
los over all 37 available snapshots from z = 2 to z = 0 
to determine the redshift evolution. We determine the 


http://www.mpa-garching.mpg.de/galform/millennium-If/ 
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number of subhalos N s h inside the search radius for each 
host halo (identified with a Friends-of-Friends ID num¬ 
ber), which can be either empty (iV s h = 0 , isolated host 
halo) or not empty (iV sh > 0, non-isolated host). 

The Millenium-II simulation is performed for a flat 
ACDM cosmology with parameter^ D m = 0.25, D& = 
0.045, D a = 0.75, h = 0.73, a 8 = 0.9, n s = 1. The 
size of the box is L = 100/i -1 Mpc (comoving) on a side. 
The spatial resolution is l/i - 1 kpc, so the search radius 
for subhalos is not affected by resolution effects. The 
mass resolution of the simulation is 6.89 x 10 6 /i _ 1 M 0 , 
so effectively our threshold on the mass of the satellite 
halo is verified by all subhalos resolved in the simulation, 
regardless of their mass. The host halo is required to have 
30,000 to 100,000 particles, which corresponds to the halo 
mass bin stated above. 

No public N-body simulations with the required mass 
resolution exist for other models than the “vanilla” 
ACDM. Fortunately, numerical simulations of one cos¬ 
mology can be re-mapped into a different one accurately 
both for the dark matter field and the corresponding ha¬ 
los [281(29] , Using this algorithm we have re-mapped the 
results of the Millennium-II ACDM simulation to cos¬ 
mologies with other values of A. Note that by keeping 
the geometry fixed (fl m + fl\ = 1) and the early universe 
quantities (such as the physical matter density D m /i 2 , the 
amplitude of primordial perturbations, and the baryon 
fraction) fixed the current value of Da specifies the cos¬ 
mology. 

It is challenging to keep track of changes of all other 
cosmological parameters and it is clear that other changes 
could mimic some of the effects that we emphasise here, 
for example significant amount of massive neutrinos 
would also suppress the growth of structure. However, 
it has been shown [30] that even modest modifications to 
a range of cosmological parameters leads to adverse con¬ 
sequences for the abundance of life. As such, we focus 
on changing only the value of the cosmological constant, 
while keeping the conditions of the early universe fixed. 

Fig.[l]shows the number of isolated halos per comoving 
1000 Mpc 3 for four relevant flat cosmologies spanning the 
range 0 < Pa/M^ < 2.7 x lO" 123 : 1) D m = 0.15, D A = 
0.85, pa/A/pj = 2.7 x 10“ 123 , 2) (fiducial ACDM) D m = 
0.25, D a = 0.75, pa/M^ = 1.4 x lO" 123 , 3) D m = 0.5, 
Da = OAjPa/A/pj = 4.5 x 10~ 124 and 4) (Einstein de 
Sitter) D m = 1, Da = 0, p\ = 0 . In the A dominated 
models, only fairly recently, on a cosmological time scale, 
the number of GRB-protected halos grew significantly. 
In fact, a large numbei[^] of isolated halos appear only 
below z « 1.5. This corresponds to a lookback-time of 


9 Unless otherwise stated, cosmological parameters are at z = 0. 

10 For a MW-like halo in the concordance ACDM model, 1/3 of the 

halos fulfil the 50 kp isolation criterion, so our MW is not a very 

special halo 
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FIG. 2. The probability distribution function for pA nor¬ 
malised to the probability for the Einstein-de-Sitter (EdS, 
A = 0) case. The solid (dashed) line corresponds to exclusion 
regions of 20 (50) kpc and result from this work. For high 
A values, the exponential suppression of Ref. m takes over. 
This is indicated by the shaded regions. The arrow shows 
the value of pA for the concordance ACDM model. Our re¬ 
sult suppresses the probability of low values of A previously 
allowed, and in fact favoured, by the argument in mm- 


about 7 Gyr, not dissimilar from the age of Earth. As 
structure formation proceeds faster in the A = 0 universe, 
the amount of sub-structure grows faster and earlier: the 
number of isolated halos is much smaller than in the case 
of a A dominated universe. 

In order for a particular galaxy to harbour life, it must 
reside within a habitable region of the parameter space. 
The halo mass should be in the range 10 11 ' 5 < M/Mq < 
10 12 as to give rise to a large galaxy like the Milky Way 
so that it has significant outer regions in which the GRB 
rate is low but there are sufficient number of stars. Ad¬ 
ditionally the galaxy should undergo sufficient chemical 
evolution so that its average metallicity is relatively large 
thereby reducing the GRB rate. Furthermore, the halo 
should be young enough that stars with mass M m j n have 
not left the main sequence. For M > 0.5Mg stars this 
corresponds to (< 50 Gyr), for M m j n > 0.7 Mq this cor¬ 
responds to < 20 Gyr. For example for the model with 
Da = 0.5 the age of the Universe is ~ 20 Gyr, by the 
time Da takes over the expansion, the universe is already 
too old and out of the habitable epoch for all main se¬ 
quence stars above M min = 0.7 Mq. For M m - m = 0.5 Mq 
we estimate this to correspond instead to Da = 0.2 
( Pa/M £ = 1.2 x 10~ 124 ). 

One can imagine waiting for a long time for the GRB 
rate to be sufficiently low. The current decrease of the 
GRB rate with time [16] at z < 1 is much flatter than the 
decrease in star formation rate (see their Fig. 9); extrap¬ 
olating this rate implies that the star formation will be 
exhausted sooner than GRBs. By z ~= -0.4- when flat 
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universes with a cosmological constant that is 0.3-0.6 of 
that in the vanilla model become similar to the present 
Universe with Ua = 0.75, the GRB rate will be down 
by a factor 2, while the star formation rate will be down 
by a factor 10. Since the decrease in star formation rate 
is the dominant effect, for simplicity in what follows we 
assume that the GRB rate remains constant. This ap¬ 
proximation does not change qualitatively our argument. 

In the range of p\ that we have explored, and for 
M m in = 0.7 Mq, we find that the number density of hab¬ 
itable and isolated halos can be approximated bjJ^j] 

I{Pa) ~ KPA + V- (!) 


The coefficients of the fit may change for different val¬ 
ues of Af m i n or different radii but the qualitative be¬ 
haviour remains with a sharp decrease for Ua < 0.2 
(/5a < 1-2 x 10 -124 ). For larger values of A this rela¬ 
tion must flatten since virtually all MW-size halos are 
isolated for > 0.85. We also note that, always in 
the range of p\ explored, the number of MW-sized halos 
in the same volume is roughly constant, implying that 
I {pa) is roughly proportional to the fraction of isolated 
halos (for larger values of A when this fraction reaches 
unity it is assumed to become a constant). If we assume 
that, for the values of p\ we considered so far, /(pa) 
is proportional to the probability p{I\p\), to infer the 
posterior probability for having a cosmological constant 
p(A|/) oc p(A)p(/|A) we must define a prior distribution 
p{ A). Previous studies take a flat prior on pa based on 
the premise that its value may arise from the cancella¬ 
tion of much larger terms. However given the large un¬ 
certainty regarding the nature of dark energy, this might 
not be the case. A more conservative approach is to se¬ 
lect an uninformative prior such that p(A) oc 1 /pa (see 
also El)- In this case one shall not consider negative val¬ 
ues of the cosmological constant, although they are likely 
to generate even more inhospitable consequences. 

The other mechanism by which A influences the habit¬ 
ability of the Universe is the suppression of galaxy forma¬ 
tion. This was studied in detail in [50; who used the frac¬ 
tion of baryons F{p) residing in halos of mass p = £ 2 M, 
where £ is the matter density per photon, as a proxy for 
the habitability of the universe. They present the follow¬ 
ing prescription for the late time solution 


F{p) = erfc 


^1 /3 

£ 4 / 3 Qs(p)Goo 


( 2 ) 


Here A and can be considered fixed quantities, A = 
5.59 and = 1.43; we take s = 28 and £ 4 Q 3 = 10 -124 


11 For 50 kpc radius r/ = 0.1 and k = 1.1 X 10 123 for pa in Mp L . 
For 20 kpc radius we find that a broken power law is a slightly 
better fit: rei = 0.72 X 10 123 , pi = 0.52; re 2 = 0.24 X 10 123 , 
r ]2 = 0.72 respectively at low and high values of p\. 


following [3D] (see e.g., their Fig 7.). It is clear from 
Eq-1 and e.g., Fig. 7 of [30] that this imposes a sharp 
suppression for pa > 10~ 122 Mpj. 

In summary, 

p(A|/)cxp(A)F(A)/(A) (3) 

where F{ A) gives the number of suitable MW-size halos 
and describes the sharp upper cut-off imposed by eq. [2] 
and /(A) imposes a much slower suppression towards 0. 
In Fig. [2] we sketch the probability distribution function 
for pa, normalised to the Einstein-de-Sitter case, using 
the computed abundance of isolated Milky Way-size ha¬ 
los from Fig. 1. The low-end suppression is described 
by Eq. 0 For high values of pa the suppression dis¬ 
cussed in (H El ED], effectively due to the requirement 
of cosmic structure to form, becomes relevant. It is not 
surprising then that our universe has a value of A « 1 
(pA ~ 10~ 123 A/p[). It is instructive to compare to Fig.7 
of f3D] and notice that the allowed region for pa is now 
greatly reduced. 

Since on average there is one isolated region in a patch 
of 10 Mpc radius (which, incidentally, is the mean in¬ 
ter galaxy cluster separation today), to ensure (at ~ 3<r) 
that there is at least one habitable galaxy in the observ¬ 
able universe, the horizon size should be at least 40 Mpc 
i.e. not more than ~ 100 times smaller than the cur¬ 
rent horizon size. Given our specifications for galactic 
habitability, both in terms of the required separations 
between galaxies, and the minimum age of the Universe 
which permits the formation of planets, a large universe 
is necessary for life to emerge. A loitering model with a 
finely tuned cosmological constant could satisfy the age 
requirement, but does not provide adequate inter-galactic 
spacing. 

If our location in the Universe, and potentially the mul- 
tiverse, is preferentially selected by the absence of nearby, 
GRB-hosting small halos, then we will expect to find that 
most Milky-Way sized galaxies beyond the Local Group 
possess a heightened number of LMC-like satellites. To 
some extent this has already been observed [3T]. The 
reduced number of observed Milky-Way satellites com¬ 
pared to predictions derived from simulations is often at¬ 
tributed to the inefficiency of star formation within low 
mass halos. Yet this puzzle may also be partially resolved 
by our proposed selection effect. 

Another interesting implication is that if the ampli¬ 
tude of fluctuations Q is increased, then halos of a given 
mass form earlier, but they do so in a very crowded en¬ 
vironment. In previous investigations Q is one of the 
few cosmological parameters which could be enlarged by 
an order of magnitude without any clearly adverse effects 
(e.g., Fig. 12 in [30]). Within the context of GRBs-limited 
habitability, there is likely to be much less freedom in this 
parameter. 

In summary, we have shown that A plays a crucial role 
at creating habitable regions for galaxies in a habitable 
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epoch. These considerations may be used to disfavour 
very low values for A. Negative values of A will yield 
even more satellites and hence these arguments strongly 
disfavour such values. 
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